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Measurements of the temperature of Mo shocked to ~60 GPa and then released to ^^28 GPa 
were previously attempted using high explosive driven flyer plates and pyrometry. Analysis of the 
radiance traces at different wavelengths indicates that the temporal evolution of the radiance can 
be explained by a contribution from the LiF window to the measured thermal radiation. 

Fitting the radiance traces with a simple model, supported by continuum dynamics studies which 
were able to relate structures in the radiance history to hydrodynamic events in the experiment, the 
contribution of the window was obtained and hence the temperature of the Mo sample. The shock- 
and release temperature obtained in the Mo was 762±40K which is consistent with calculations 
taking the contribution of plastic work to the heating into account. The radiance obtained for the 
LiF window shows a non thermal distribution which can be described by a bulk temperature of 
624± 112K and hot spots (less than 0.5% in total volume) within the window at a temperature of 
about 2000K. 

PACS numbers: 62.50.Ef, 07.20.Ka 
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I. INTRODUCTION 

The measurement of temperature in shock physics ex- 
periments is of paramount importance since temperature 
is required as constraint (additional to pressure and vol- 
ume) for the development and credibility of robust multi- 
phase equation of states (EOS). The application of mod- 
ern material strength models requires these EOS with 
accurate temperatures and phase boundaries, particu- 
larly melt curves. Temperature is an important parame- 
ter in geophysics and planetary sciences for an accurate 
description of planetary structure and astrophysical im- 
pacts. Many measurements of the EOS of planetary ma- 
terials, and the pressure calibration on which diamond 
anvil cell studies rely, are made using shock experiments, 
and rely on temperature corrections to infer states off 
the shock Hugoniot. General research and model devel- 
opment for the response of materials to extreme dynamic 
conditions is focusing on multi-scale approaches, in which 
physics-based models (as opposed to empirical relations) 
for thermally-activated processes including plastic flow, 
phase changes, and chemical reactions requires an ac- 
curate knowledge of the temperature. As a specific ex- 
ample, the melting curve and kinetics of Beryllium and 
Carbon (Diamond) at elevated pressures are important 
to understand the behavior of ablator capsules for the 
national ignition facility (NIF) to succeed in harvesting 
fusion energy from laser driven capsule implosions^. 

Pyrometry is a form of thermal emission spectrom- 
etry in which the emission spectrum is collected in a 
small number of spectral bands of significant bandwidth, 
compared with energy-dispersive spectrometers such as 
prisms and gratings. Energy-dispersive spectrometers 



are not practical for most shock experiments at temper- 
atures below a few thousand Kelvin because the thermal 
emission is too weak to take advantage of the relatively 
fine spectral resolution. Pyrometry is thus the most 
promising technique applicable to many materials and 
experimental configurations (besides Neutron Resonance 
Spectroscopy^!^ and Raman Spectroscopy'', both with 
their advantages and limitations) to achieve the goals 
stated above. Nevertheless, although pyrometry has been 
fielded on dynamic loading experiments for more than 
four decades''' it still suffers from problems such as back- 
ground light, thermal and non-thermal emission from 
the window material and sample/window interface effects 
which are hard to take into account. 

Neutron Resonance Spectroscopy (NRS) has been in- 
vestigated as a possible bulk-temperature diagnostic 
for shock physics experiments, capable of measuring 
the temperature within opaque samples. Initial ex- 
periments on shock loaded Mo showed significant dis- 
crepancy with theoretical predictions. To investigate 
this discrepancy, microsecond-duration pyrometric mea- 
surements were performed on Mo samples through LiF 
windows. These experiments used a shock generation 
method that was the same as used in the NRS experi- 
ments, and also included Doppler velocimetry measure- 
ments of the surface velocity history of the sample to 
verify the loading conditions. Previous analysis of the 
pyrometry experiments^ did not succeed in extracting 
temperatures. Here we report a more rigorous study 
of the effects of shock and release waves on pyrometry 
measurements with a release window, identifying two- 
temperature population from pyrometry data. We are 
able to correlate features of the radiance history at a level 



usually dismissed as too complicated for further analy- 
sis and provide valuable insight into dynamic processes 
occurring in shock experiments using a window. This 
allows us to extract temperatures more accurately and 
with greater confidence. The results presented here are 
the first simultaneous extraction of sample and window 
temperature to our knowledge which enhances the ability 
to reliably interpret pyrometry data. This is a first step 
and key to future research and physics-based understand- 
ing of response of condensed matter subject to dynamic 
loading and heating. 



II. EXPERIMENTAL SETUP 

As in the NRS experiments, a shock was induced in 
the Mo sample by the impact of an Al disk accelerated 
to a speed of about 3.5 km/s by means of high explosive 
(HE) gases-. A LiF window was glued to the backside of 
the Mo disk (where the measurement was made) in or- 
der to maintain high pressure over an area accessible for 
contactless temperature measurements^. Emitted ther- 
mal radiation was focused onto a 1-mm core-diameter 
near-infrared (NIR) fiber as well as onto a 1-mm core- 
diameter visible-glass, low OH fiber. Two different multi- 
channel pyrometers were used to infer thermal radiance 
over a wide range of wavelengths: 5 channels in the vis- 
ible and NIR wavelength region, using photomultiplier 
tubes (PMTs), and 4 channels in the NIR, using InSb 
detectors. The velocity history at the Mo/LiF interface 
was measured by laser Doppler velocimetry, using a Ve- 
locity Interferometer System for Any Reflector (VISAR)^ 
in order to infer the pressure history applied to the sam- 
ple. 

In the case of low-temperature pyrometry measure- 
ments, special precautions have to be taken in order to 
avoid background light, which can easily overwhelm espe- 
cially the short wavelength channels. These backgrounds 
can be generated either by bright HE gases blowing by 
the target, air lighting up due to shock luminescence from 
nitrogen, ejecta or jets generated due to improper surface 
preparation, or sharp edges. In order to minimize the 
background, the experiments were performed in a vac- 
uum (10-3 Torr). Two CaF lenses were used to focus 
the thermally emitted light onto a 1-mm diameter NIR 
fiber^^ (centered 2.5 mm off axis at the fiber bundle) and 
a 1-mm diameter visible low OH fiber— (centered 2.5 mm 
diametrically opposed). The center fiber was a VISAR 
probe^^ used to determine the velocity of the sample af- 
ter impact. The LiF window, 30 mm diameter by 20 mm 
thick, was included to maintain an elevated pressure in 
the Mo when the shock reached its surface, avoiding a 
release to atmospheric pressure. The LiF was attached 
to the Mo coupon using Loctite 326 glue, which has been 
previously found^" not to cause considerable amounts of 
thermal emission. This simplifies relating the surface 
temperature to the bulk temperature of the shocked sam- 
ple. However thermal light emission from both the glue 
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FIG. 1: Schematic setup of the experiment. S, molybdenum 
sample (064 x 6 mm); W, LiF window (025 x 20 mm), F; 
aluminum flyer (064 x 5 mm, 20 mm gap between F and S); 
HE, 9501 high explosive (7 mm gap between HE and F); L, 
CaF2 lenses (025 mm, f=50 mm); OF, optical fibers (01 mm 
C2 fiber, 01 mm low OH fiber & VISAR fiber bundle). 

layer and the LiF is a concern. 

The NIR pyrometer collimates and divides the incom- 
ing light from a single NIR fiber into four spectral ranges, 
using three custom dichroic beamsplitters. These four 
coUimated beams are spectrally narrowed by bandpass 
filters centered at 1.8 ^m, 2.3 /xm, 3.5 /j,m and 4.8 /xm and 
then refocused onto the 1 mm2 active area of the 50 MHz 
bandwidth InSb detectors using ZnSe lenses. The lower 
temperature limit of the two longer wavelength channels 
is about 340 K (assuming an emittance of 1). For more 
information on the NIR instrument and on data analysis 
see Ref. [nl 

The 5 channel visible/NIR pyrometer also uses dichroic 
beamsplitters to spectrally divide the incoming light into 
five beams, which are then refocused onto the active area 
of photo multiplier tubes (PMT). Also using bandpass 
filters, these five channels are centered at 505 nm, 725 
nm, 850 nm, 1.23 /xm and 1.59 /zm. A holographic notch 
filter is used to suppress the bright 532 nm laser light 
used for the VISAR measurements. More information 
on this instrument can be found in Ref. [l^ . 

For further details on the HE driven flyer system see 
references 0and[il. 

III. RESULTS 

Eight experiments have been performed at the 
"Chamber-8" high explosive experimental facility at the 
Los Alamos National Laboratory in September 2004; 
these experiments are summarized in Table HI 

Figure [2] shows the radiance as function of time as well 
as the sample/window interface velocity for experiment 
^06. The time of shock breakout is set to t = fis. If 
no unwanted background light occurs the radiances are 
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TABLE L Experiments performed. Surface condition 
(free surface or window), diagnostics (V, VISAR; IR, IR- 
pyrometer; vis, visible-pyrometer) 



Exp# 


Surface 


Diagnostics 


Comments 


01 


LiF window 


V, IR 


Data clipped 


02 


LiF window 


V, IR 


Data lost 


03 


LiF window 


V, IR 


Good data 


04 


free surface 


V, IR 


HE problems, no data 


05 


free surface 


V, IR 


Good data 


06 


LiF window 


V, IR, vis 


Good data 


07 


LiF window 


V, IR, vis 


Good data 


08 


free surface 


V, IR, vis 


Good data 
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FIG. 2: Measured radiances for seven wavelengths and sam- 
ple/window interface velocity for experiment #06. 



expected to be constant from the time of shock breakout 
until the release wave from the back of the flyer reaches 
the sample/window interface (about 0.7 /iS after break- 
out). At this time the sample/window interface is de- 
celerating to a lower velocity as can be seen from the 
particle velocity trace. 

Continuum dynamics simulations were used to inves- 
tigate the origin of the varying signals in the shocked 
state by comparing them with the various hydrodynamic 
events occurring in the experiments. One and two di- 
mensional (ID and 2D) simulations were performed, us- 
ing general- purpose multi-physics hydrocodes. Shock dy- 
namics calculations were also performed using the same 
material properties, solving the Rankine-Hugoniot equa- 
tions to obtain shock states and hence wave speeds more 
precisely than by derivation from the discretized hy- 
drocode solutions^. The ID and 2D hydrocodes both 
used a finite difference representation of spatial fields, a 
second-order predictor-corrector algorithm for time inte- 
gration, and artificial viscosity to stabilize shocks. The 
ID simulations, along the axis of symmetry, were La- 
grangian, avoiding inaccuracies caused by numerical ad- 
vection, and used the LAGCID hydrocode— . The 2D 
simulations, axisymmetric in the axial-radial plane, were 
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FIG. 3: x-t-diagram (lower part: black lines, boundaries be- 
tween flyer (F), sample (S) and window (W); grey lines, shock 
front; dashed lines, head and tail of flyer release waves), and 
measured spectral radiance at 2.3 fim for experiment 06 (up- 
per part: the time has been shifted so that the breakout at 
the sample-window interface occurs at the same time as at 
the simulated x-t-diagram). 



Eulerian for robustness in treating highly distorted flow 
from sharp corners, with an operator-split third order 
van Leer flux limited method for advection, and used the 
EUL2D hydrocodei^. 

The simulations included the HE-driven flyer, the Mo 
target, and the LiF window. The HE acceleration sys- 
tem itself was not modeled: the flyer was treated as flat 
and moving at a constant speed, starting at the instant 
of impact. In reality, the flyers were slightly dished, still 
accelerating slightly, and reverberating, but these details 
should not make a signiflcant difference to the pyrom- 
etry data. The equations of state (EOS) used the cu- 
bic Grueneisen form with published parameters for each 
material^^. Some simulations were repeated using EOS 
from the SESAME tabular compendiumA*.^; the difference 
was negligible. Strength was treated using the Steinberg- 
Guinan modeU^, with published parameters for each 
material^^. In the ID simulations, the cell size in the 
flyer, sample, and window was 0.1, 0.05, and 0.2 mm re- 
spectively. In the 2D simulations, the cell size was 0.1 
mm. 

A position-time diagram for the key mechanical waves 
from the simulations is shown in figure [3] (lower part) 
together with the measured spectral radiance at 2.3 /im 
for experiment #06. 

The radiance increased from the time of shock break- 
out until about 2.1 fis. Then the radiance decreased at a 
rate slightly lower than the increase. The kink in the ra- 
diance at about 3.1 /is occurs at about the time when the 
release from the edge of the window reaches the center 
of the sample-window interface, cooling it by adiabatic 
expansion. The shock wave reaches the rear surface of 
the window at about 3.75 fis, at this time the signal of 
all wavelengths goes into saturation because of a bright 
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FIG. 4: Schematic of the window-sample combination. S, 
sample; Ws, shocked window; Wu, unshocked window, R1-3, 
interface reflectivities. 



flash that occurs at this event. 

The recorded pyrometer output signals for experiment 
#01 have been clipped because the settings of the dig- 
itizer have been too sensitive (the increase of radiance 
with time was not expected), the measured radiances for 
the other successful experiments look very similar to the 
ones shown in figure [2] (for experiment #06) and are not 
shown here. 

The radiance traces for the free surface experiments are 
not analyzed and discussed here because, additionally to 
problems with background light, the surface temperature 
was not homogeneous. (This can be inferred from the re- 
sponse of different wavelengths shortly after the time of 
shock breakout. For more details on spatial temperature 
non-uniformities see reference [20I) The background light 
at the free surface experiments was most likely caused 
by thermal radiation from APIEZON Q©, a soft, black, 
putty-like substance that was intended to absorb ther- 
mal radiation from hot jets formed at the edge of the 
sample/window interface. 



IV. DATA ANALYSIS 

Figure U] shows the parts involved in the process of 
emission of thermal light for the sample/window com- 
bination used in these experiments, including the effect 
of the interface between shocked and unshocked window 
material. 

Thermal emitted radiance from the sample (Lgampie) 
is multiply reflected between the interfaces Ri, R2 and 
R3 . These reflectivities can be obtained from the optical 
properties of the sample, the shocked and the unshocked 
window. The index of refraction of the shocked window 
can be estimated from the Gladstone-Dale relational, 
from this it can be seen that R2 is very small and can be 
neglected'^'^ . The index of refraction and extinction coef- 



flcient of the unshocked window can be found in Ref. |22, 
the absorption in this part can be neglected as well. The 
shocked part of the window can emit thermal radiation 
(J-'window) as well as absorb thermal radiation emitted by 
the sample. The emissive power of translucent materi- 
als (equivalent to emissivity of opaque materials) can be 
expressed by (1 — e~"^)^, and the absorption can be ex- 
pressed by e~"^ , where x is the thickness of the shocked 
layer and a is the absorption coefficient. The absorption 
coefficient is related to the extinction coefficient by the 
following relation: 



A-TT ■ k 

An 



(1) 



where k is the extinction coefficient and Aq is the wave- 
length in vacuum. Thus the radiance from this assembly 
can be expressed as: 



L = 



(1 - R3) 



sample ' ^ ' 



(1 - i?i • i?3 • e-2"^) 

+ L^^ndo^o ' (1 + i?l ' e""^) • (1 - e""^)] (2) 



where e is the effective emissivity of the sample into the 
window, taking the interface effects and absorption of the 
glue layer into account. The emissivity is mainly deter- 
mined by the optical properties (n and k) of the sample 
and the window, but can increase if the interface rough- 
ness increases following shock breakout, for instance if 
the interface is Richtmyer-Meshkov unstable or because 
of flow or jetting from surface texture. 

The first term in equation^] ( 1 - i?3 ) / ( 1 - i?i • i?3 • e-^"^ ) 
is very close to 1 (between 0.98 and 1.0 in all practical 
cases) and can be set to 1 without compromising the ac- 
curacy of the analysis. The term (l-|-i?i-e~"^)-(l — e~"^) 



can be simplified to (1 



'^) without compromising 



the fitting results for the sample radiance, introducing an 
error in the window radiation and absorption coefficient 
of less than 10%, which is far below the uncertainty of 
the fitting parameters. Thus equation [2] can be approxi- 
mated as: 



sample ' * 



~^ ^window ' (1 6 ) (3) 



The thickness x of the shocked layer is given by: x — 
{us — Up) ■ t , where Ug is the shock velocity in LiF (7.10 
mm//is at a pressure of 27GPa) and Up is the particle 
velocity at the sample/window interface (1.45 mm//is on 
release into LiF from a Hugoniot pressure of 64GPa in 
Molybdenum) and t is time after shock breakout at the 
Mo/LiF interface. Equation[3]was used to fit the radiance 
traces of experiments #01, 03, 06 and 07. The results are 
shown in table |TT1 

While the emissive power of the window material'^ was 
obtained from the fitting parameters, and hence the tem- 
perature of the shocked window can be obtained directly, 
the inferred temperature of the sample depends strongly 
on the effective emissivity e of the sample-window com- 
bination. Since this property is very hard to determine 
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TABLE II; Fitting results for the experiments #01, 03, 06 
and 07. Where there are no data available (due to low signal 
level or poor quality of the fit) the cells have been left blank. 



Exp# 


A 


e • Lmo 


LhiF 


a,LiF 




(ptm) 


(W/m^ ■ nm ■ sr) 


{W/m? ■ nm ■ sr) 




01 


3.5 


0.117 


2.27 


0.0109 


4.8 


0.051 


0.229 


0.0988 


03 


3.5 


0.110 






4.8 


0.062 








0.85 






0.38694 




1.23 




16.3 


0.11307 




1.59 




15.8 


0.08272 


06 


1.8 


0.072 


4.73273 


0.0472 




2.3 


0.088 


3.2168 


0.0239 




3.5 


0.06401 


0.85708 


0.03602 




4.8 


0.02808 


0.45 


0.11792 




1.23 


0.004 








1.59 


0.025 






07 


1.8 


0.102 






2.3 


0.108 








3.5 


0.132 








4.8 


0.069 




0.1234 



experimentally, it is common practice to assume lower 
and upper emissivity bounds and use these bounds to 
calculate upper and lower bounds for the temperature^!. 
This has been done for experiments #06 and 07 (where 
a sample radiance T-isampie was obtained) with the fol- 
lowing lower 6/ and upper e„ emissivity bounds: A=:1.23 
/im: ei = 0.2, e„ = 1.0; A=1.59 /xm: e/ = 0.15, e„ — 
0.8; A=1.80 /im: ei = 0.12, e„ = 0.8; A=2.30 /im: = 
0.1, e„ = 0.6; A=3.50 iim: ei = 0.08, e„ = 0.4; A=4.80 
lim: ei = 0.05,6^= 0.3. These estimates are based on 
the room temperature values (obtained from the opti- 
cal properties) for polished surfaces^ii^ and the melting 
temperature values at ambient pressures^!. The calcu- 
lated true temperature bounds for experiment #06 and 
07 as function of wavelength are shown in figure O (for 
readability reasons the wavelengths for experiment #06 
have been shifted by 0.1 to the red). 

It can be seen that the uncertainty due to unknown 
emissivity is far higher for longer wavel eng ths, as ex- 
pected. (For an explanation see reference l20l) The tem- 
perature ranges for all wavelengths overlap between 722 
K and 802 K; this is the range where the temperature of 
the sample shocked to a Hugoniot pressure of 63.9±2.4 
GPa and released into LiF (at an interface pressure of 
27.1±1.0 GPa) is determined. The Hugoniot pressures in 
the sample have been obtained from the measured parti- 
cle velocity at the sample window interface, the pressures 
in the LiF window have been obtained from jump condi- 
tions and known shock properties of the window material, 
through the published EOS and strength modelii. 
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FIG. 5: Upper and lower temperature bounds as function 
of wavelength for experiments #06 and 07 (the wavelengths 
for experiment #06 have been shifted 0.1 nm to the red for 
readability reasons). 
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FIG. 6: Apparent window temperature as function of wave- 
length (squares) from experiment #06. The solid line is the 
result from fitting the measured data with a simple two tem- 
perature mode&, the results of the fitting process are pre- 
sented in the text below. 



The temperature of the LiF window as function of 
wavelength from experiment #06 is shown in figure [HI 

If the window is at a homogenous temperature and no 
background light occurs, the measured window tempera- 
ture should be independent of wavelength. The increase 
of the measured window temperature (apparent tem- 
perature) with decreasing wavelength indicates a non- 
homogenous window temperature. Fitting the results 
with a simple two temperature model (for details see ref- 
erence[20|) gives a temperature of the window of 624±112 
K with hot cells at a temperature of 1945±210 K. The 
area fraction (as seen from the pyrometer) of these hot 
cells was determined to be 3.0-1-1.5%, corresponding to a 
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FIG. 7: Absorption coefficient for shocked LiF (27.1±1.0 
GPa) as function of wavelength for aU wavelengths and ex- 
periments where these data could be obtained from the fit- 
ting process. The solid line shows literature values for LiF at 
standard conditions. 



volume fraction of less than 0.5%. The absorption coeffi- 
cient as function of wavelength for the shocked LiF win- 
dow is shown in figure [T] The solid line shows data from 
reference [22I at room temperature and ambient pressure. 
It can be seen that the absorption coefhcient increases by 
roughly two orders of magnitude compared to the litera- 
ture data independent of wavelength. 

It can be seen from figure O that the radiance increases 
as long as the thickness of the shocked layer in LiF is in- 
creasing (from shock breakout at the Mo-LiF interface 
at 1.0 until the release wave from the back of the 
aluminum flyer reaches this interface at about 2.1 /is). 
The thickness of the shocked LiF layer decreases after 
t=2.1 ^s because the release wave propagates at a higher 
velocity than the shock wave and hence the measured ra- 
diation decreases between t=2.1 /is to about 3.7 ^s. At 
3.7 /iS the shock front in the LiF reaches the free sur- 
face of the window, and the measured radiance increases 
abruptly, which is assumed to be due to non thermal light 
emission when the window fractures on release from the 
shocked state and possibly the air shock. 

From fitting the spectral radiance traces at different 
wavelengths and assuming lower and upper bounds for 
the effective emissivity of the sample/ window assembly, 
a sample temperature (shocked to 63.9±2.4 GPa and re- 
leased to a pressure of 27.1±1.0 GPa into LiF) at the 
interface of 762±40 K could be obtained. The release 
temperature into LiF was calculated to be 670±25 K us- 
ing the Steinberg-Guinan strength model to predict the 
contribution to heating from plastic work [25|. The re- 
sulting temperatures are in relatively good agreement^ 
with experiments performed on a gas gun at a slightly 
lower pressure (T=683±41 K at a Hugoniot pressure of 
58.7 GPa, with a release pressure of 24.8 GPa into LiF). 
A summary of all experimental and calculated results can 



FIG. 8: Measured and calculated temperatures for molyb- 
denum releasing into LiF for two different release pressures. 
Open square, this experiment; full square, gas gun experiment 
[2a |: circles, calculations taking strength into account via the 
Steinberg-Guinan model [l^; triangles, calculations ignoring 
strength effects [l9|. 

be seen in figure [H 

The measured temperature in the LiF window is 
624±112 K with hot spots (3.0±1.5% in area) with a 
temperature of 1945±210 K. These hot spots are be- 
lieved to occur due to the slightly non uniaxial loading 
of the sample and window [3,] due to the dishing of the 
flyer plate which most likely causes local heating by shear 
strain localization. This localized heating in LiF windows 
was also observed in infrared imaging experiments [26| . 
The temperature of the LiF window was calculated to be 
569±15 K [25II which is in excellent agreement with the 
temperature inferred from the experiment. An analysis of 
the free surface experiments was not possible because of 
temperature non uniformities combined with background 
light caused by the hot APIEZON Q(c)(used to block 
thermal emission from hot jets from the edges, as dis- 
cussed above). 



V. CONCLUSIONS 

Reliable pyrometric temperatures were obtained at 
shocked and released surfaces despite an o(100%) back- 
ground from shocked window material. This improve- 
ment was achieved by using the emission spectrum to in- 
fer a volume fraction of hotspots in the window, and was 
supported by relating changes in the emitted radiance to 
predictions of hydrodynamic events in the window. The 
measured temperature for Mo at a Hugoniot pressure of 
63.9±2.4 GPa released into LiF to an interface pressure 
of 27.1±1.0 GPa of 762±40 K is in good agreement with 
temperatures measured using a powder gun at a slightly 
lower pressure (58.7 GPa released to 24.8 GPa) of 683±41 
K. The measured temperatures are slightly higher than 
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calculations using the Steinberg-Guinan model [19| to 
take strength effects into account. It is not clear whether 
the remaining discrepancy is caused by additional sources 
of thermal emission not taken into account, or inaccuracy 
in the models used in the simulations. However, this level 
of agreement is unusually good for thermal emission from 
shocked metals. 

Besides the release temperature of Mo, the shock tem- 
perature of LiF at a pressure of 27.1±1.0 GPa was also 
obtained and is in good agreement with calculations. The 
shock temperature in the window was much less sensitive 
to the strength model. 



Because of problems controlling the thermal back- 
ground, free surface temperatures could not be extracted 
from the experiments without a window. 
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